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(57) ABSTRACT 

Methods for using modified single wall carbon nanotubes 
(“SWCNTs”) to detect presence and/or concentration of a gas 
component, such as a halogen (e.g., Cl 2 ), hydrogen halides 
(e.g., HC1), a hydrocarbon (e.g., C„H 2w+2 ), an alcohol, an 
aldehyde or a ketone, to which an unmodified SWCNT is 
substantially non-reactive. In a first embodiment, a connected 
network of SWCNTs is coated with a selected polymer, such 
as chlorosulfonated polyethylene, hydroxypropyl cellulose, 
polystyrene and/or polyvinyl alcohol, and change in an elec- 
trical parameter or response value (e.g., conductance, current, 
voltage difference or resistance) of the coated versus 
uncoated SWCNT networks is analyzed. In a second embodi- 
ment, the network is doped with a transition element, such as 
Pd, Pt, Rh, Ir, Ru, Os and/or Au, and change in an electrical 
parameter value is again analyzed. The parameter change 
value depends monotonically, not necessarily linearly, upon 
concentration of the gas component. Two general algorithms 
are presented for estimating concentration value(s), or upper 
or lower concentration bounds on such values, from measured 
differences of response values. 

10 Claims, 22 Drawing Sheets 
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COATED OR DOPED CARBON NANOTUBE 
NETWORK SENSORS AS AFFECTED BY 
ENVIRONMENTAL PARAMETERS 

ORIGIN OF THE INVENTION 

This invention is a Continuation-In-Part of U.S. Ser. No. 
1 1/178,079, filed 8 Jul. 2005, and was made, in part or whole, 
by one or more employees of the U.S. government. 

The U.S. government has the right to make, use and/or sell 
the invention described herein without payment of compen- 
sation, including but not limited to payment of royalties. 

FIELD OF THE INVENTION 

This invention relates to use of carbon nanotube networks 
as sensors of chemical substances. 

BACKGROUND OF THE INVENTION 

Chemical sensors have been developed for decades now to 
detect gases and vapors at various concentration levels for 
deployment in a wide range of industrial applications. The 
detection usually centers on a change of a particular property 
or status of the sensing material (such as temperature, elec- 
trical, optical characteristics, etc.) upon exposure to the 
chemical species of interest. The selection of sensing material 
itself has spanned across the periodic table with a range of 
inorganic, semiconducting elements and organic compounds 
either in bulk or in thin film form. Perhaps the most widely 
investigated class of sensors is the high temperature metal 
oxide sensor, due to its high sensitivity with tin oxide as an 
example of sensor material. The most common Sn0 2 sensor 
platform is a chemiresistor wherein the transport character- 
istics of a conducting channel of tin oxide are modulated by 
the adsorption of chemical species at elevated temperatures 
(Ti^350° C.). Other types of sensors include electrochemical 
cells, conducting polymer sensors, surface acoustic wave sen- 
sors and catalytic bead sensors. 

While commercial sensors based on some of the above 
approaches are available, research continues on new sensing 
materials and transducer platforms for improved perfor- 
mance. Desirable attributes of next generation sensors 
include high sensitivity, in the parts per million (ppm) to parts 
per billion (ppb) range, low power consumption, room tem- 
perature operation, rapid response time, high selectivity and 
long term stability. Sensors based on the emerging nanotech- 
nology promise to provide improved performance on all of 
the above metrics compared to the current micro and macro 
sensors. Nanomaterials exhibit small size, light weight, very 
high surface to volume ratio, increased chemical reactivity 
compared to bulk materials, and mechanical stability so that a 
sensing material can be refreshed or regenerated many times. 
All these properties are ideal for developing extremely sen- 
sitive chemical sensors. 

Among the numerous nanomaterials, carbon nanotubes 
(CNTs) have received significant attention due to their unique 
electronic and extraordinary mechanical properties. Single- 
wall carbon nanotubes (SWCNTs) have an enormous surface 
area, as high as about 1600 m 2 /gm, which leads to an 
increased adsorptive capacity for gases and vapors. With all 
the atoms on the surface, SWCNTs are expected to exhibit a 
change in properties sensitively upon exposure to the envi- 
ronment. Indeed, electrical conductivity of SWCNTs has 
been shown to change reproducibly in the presence of gases 
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such as N0 2 and NH 3 . This revelation has resulted in the 
fabrication of SWCNT-based chemical sensors by several 
groups. 

The principal platform for such sensors has been a nano- 
5 tube field effect transistor (“CNT-FET”) with a single 
SWCNT serving as the conducting channel. This platfonn 
faces some serious difficulties for commercialization. First, 
the CNT-FET requires semiconducting SWCNTs for its 
operation, and selective growth of metallic versus semicon- 
io ducting nanotubes is not possible today. Second, if an in situ 
chemical vapor deposition (“CVD”) process is used in the 
device fabrication sequence, it is hard to make a single 
SWCNT grow horizontally in order to bridge a given distance 
between the source and the drain. Alternatively, one is forced 
15 to ‘pick and place’ a semiconducting SWNT from bulk 
samples. Finally, the chemical sensor market is too cost sen- 
sitive to rely on complex steps involved in CNT-FET fabri- 
cation resulting in low sensor yield and poor reproducibility. 

What is needed is an approach using suitably modified 
20 nanomaterials, such as SWCNTs, that can detect presence of 
certain gas components whose presence cannot be detected 
by any simple means. Preferably, the method should provide 
high sensitivity (detection of parts per million or parts per 
billion of the target gas), high selectivity, room temperature 
25 operation, low power consumption, high throughput and low 
cost. Preferably, the method(s) should extend to detection of 
other gas components with at most modest changes in the 
nanomaterial modification procedures. Preferably, this 
method should allow estimation of the effects of change of 
30 local temperature, change of local relative humidity and tem- 
poral drift (change of baseline and sensitivity with elapsed 
time). 

SUMMARY OF THE INVENTION 

35 

These needs are met by the invention, which provides a 
chemical sensor or sensor array for detecting presence, at or 
near room temperature, of one or more of N target gas com- 
ponents or molecules (NiH) in a gas mixture contained in a 
40 chamber, by any ambient being considered. The sensor con- 
tains a network of SWCNTs that is connected to a controlla- 
bly variable voltage difference or current source. The cham- 
ber may be closed, isolated and static; or, preferably, may 
allow gas flow-through and thus not be wholly isolated from 
45 the external environment. Alternatively, the chamber may be 
part or all of the external environment. 

In a first embodiment, the SWCNTs in the network are 
partly or wholly coated with a selected polymer, such as 
chlorosulfonated polyethylene, hydroxypropyl cellulose, 
50 polystyrene or poly vinylalcohol, and the target molecule may 
be a hydrocarbon C m H„, (e.g., CH 4 or C 3 H 7 or C 2 H 2 ), an 
alcohol C m H„,OH, a ketone (e.g., CH 3 (CO)CH 3 ) or an alde- 
hyde (e.g., C 2 H 5 (COH). An algorithm, applicable to any 
embodiment, provides an estimate of target molecule concen- 
55 tration. 

In a second embodiment, the SWCNTs in the network are 
doped with a selected transition element, such as Pd, Pt, Rh, 
Ir, Ru, Os and/or Au. In either embodiment, a value of an 
electrical parameter, such as conductivity, resistivity, electri- 
60 cal current or voltage difference, is measured and compared 
with the parameter value for the network with no coating and 
no doping. SWCNTs are positioned between the electrodes 
using either a solution casting process in the form of a net- 
work or in situ growth using chemical vapor deposition 
65 (CVD) techniques. Polymer coating or transition element 
doping of SWCNTs allows selective sensing of certain gases, 
as demonstrated here for chlorine (Cl 2 ), HC1, CH 4 and CO x 
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vapor. The CNT sensors may be formed on a substrate, such 
as silicon, ceramic, glass and selected polymers. The sensor 
fabrication process is scalable for manufacturing products 
that include wafer scale interdigitated electrode fabrication 
and inkjet deposition of SWCNTs and polymers for coatings 
and metal nanoparticles for doping. Coating and doping may 
be collectively referred to as “loading.” 

In a third embodiment, an array of substantially identical 
SWCNTs are (i) coated, as in the first embodiment, or (ii) 
doped, as in the second embodiment, a temperature gradient 
is imposed on the array, and a value of an electrical parameter 
is determined for each coating or doping tested, and a change 
in sensitivity is determined or estimated as a function of the 
local temperature value or gradient. 

In a fourth embodiment, a sequence of substantially iden- 
tical SWCNTs are (i) coated, as in the first embodiment, or (ii) 
doped, as in the second embodiment, different relative 
humidity values (e.g., 0, 15, 30, 50, 70 and 90 percent) or 
other environmental parameter values are imposed on each 
SWCNT sensor, a value of an electrical parameter is deter- 
mined for each coating or doping tested, and a change in 
sensitivity is determined or estimated as a function of changes 
in the environmental parameter value. 

In a fifth embodiment, an array of substantially identical 
SWCNTs are (i) coated, as in the first embodiment, or (ii) 
doped, as in the second embodiment, a value of an electrical 
parameter is determined for a selected local temperature and 
a selected local relative humidity, at the end of each of a 
sequence of selected time intervals, with temporal lengths 
from a few seconds to six months or more, if desired, and a 
change in baseline and sensitivity or “drift” is determined or 
estimated as a function of elapsed time. An algorithm is 
presented that can compensate for this drift in baseline and 
sensitivity, as a function of time. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates use of an array of sensors, each using a 
network of SWCNTs to connect two electrodes according to 
the invention. 

FIGS. 2 and 3 graphically illustrate sensor response (elec- 
trical current) versus concentration (C) for Cl 2 detection and 
for HC1 detection, using SWCNTs coated with chlorosul- 
fonated polyethylene and hydroxypropyl cellulose, respec- 
tively. 

FIG. 4 compares relative response of uncoated SWCNTs 
for presence of N0 2 , NH 3 , benzene, nitrotoluene, Cl 2 and 
HC1. 

FIG. 5 compares relative response of SWCNTs coated with 
chlorosulfonated polyethylene or with hydroxypropyl cellu- 
lose for presence of N0 2 , Cl 2 and HC1. 

FIG. 6 is a flow chart for practicing the invention using one 
or more selected coatings (e.g., chlorosulfonated polyethyl- 
ene or hydroxypropyl cellulose). 

FIGS. 7 and 8 are flow charts of algorithms for estimation 
of concentrations of gas components and for determination of 
presence of a particular gas component, where a coated or 
doped SWCNT network is used. 

FIG. 9 graphically illustrates recovery time of a measured 
electrical parameter (here, current) to presence of ultraviolet 
light. 

FIG. 10 graphically illustrates variation of a measured 
electrical parameter with time for different concentrations. 

*FIGS. 11 and 12 graphically illustrate behavior of a mea- 
sured electrical parameter and response time with varying 
temperature T. 
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*FIG. 13 graphically illustrates behavior of a measured 
electrical parameter with varying relative humidity RH, for a 
SWCNT coated with hydroxypropyl cellulose. 

*FIG. 14 graphically illustrates variation of an approxima- 
5 tion function for a measured electrical parameter with varying 
relative humidity RH for an SWCNT coated with an unspeci- 
fied substance. 

FIGS. 15-1 through 15-3 graphically present sequences of 
measurements of a response parameter W (electrical resis- 
io tance) for an array of coated SWCNTs, as a function of N0 2 
gas concentration for varying relative humidity values of 0, 
30, 50, 70, 90 percent. 

FIGS. 16A-16F graphically present measured values of a 
parameter W (electrical resistance) for a coated SWCNT, for 
15 varying relative humidity values RH and for time varying 
concentration values ofN0 2 gas. 

FIGS. 17 and 18 are flow charts of procedures for practic- 
ing the invention in the presence of temporal drift and for a 
selected environmental parameter. 

20 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

FIG. 1 illustrates a sensor 11, constructed according to an 
25 embodiment of the invention, including 32 interdigitated 
electrode fingers (“IDEFs”), 12A-1 through 12D-8, con- 
nected to first and second voltage sources, with a controllable 
voltage difference. In FIG. 1, a 1 cmxl cm chip with Q 
sensing elements and SWCNTs bridging the gaps between 
30 adjacent IDEFs (Q=2-32). Each IDE was fabricated using 
conventional photolithography with a nominal finger width of 
10 pm and gap sizes d=8, 12, 25 and 50 pm , or any other 
desired set of gap sizes. The fingers and contact pads were 
thermally evaporated Ti (20 nm thickness) and Pt (200 nm 
35 thickness) on a layer of Si0 2 , thermally grown on a silicon 
substrate. The sensing material was bulk-produced SWCNTs 
from a HiPCo process (Rice University.), purified to remove 
amorphous carbon and metal impurities according to a pro- 
cedure described in the literature. The purified SWCNTs 
40 were dispersed in a selected solvent, such as dimethyl forrna- 
mide, to create a suspension of SWCNTs. The suspension 
was sonicated, then drop -deposited onto the interdigitated 
area of the electrodes. After the solvent evaporated, the 
SWCNTs formed a network connecting two adjacent IDEFs 
45 with a corresponding voltage difference. Any residue of the 
solvent was thoroughly removed by drying the sensor under 
vacuum. The SWCNT density in the network can be varied by 
varying the amount of suspension or the concentration of 
SWCNTs in the suspension placed on the sensor area. This 
50 process produces a statistically meaningful number of 
SWCNTs bridging the gap between two terminals to give 
reproducible performance. Alternatively, the SWCNTs can 
be grown directly on the electrodes using chemical vapor 
deposition. A selected electrical parameter response value, 
55 such as conductance or resistivity or electrical current or 
voltage difference, is provided or measured for the uncoated 
and coated SWCNT networks. Differences between the sens- 
ing elements are compensated for by varying a base resistance 
located in series with each sensing element. 

60 SWCNT s yield different signal responses when expo sed to 

different gases and vapors and one must use pattern recogni- 
tion or intelligent signal processing techniques for the iden- 
tification of the gas constituent of interest. SWCNTs do not 
respond to exposure to certain gases and vapors, and in those 
65 cases, coating or doping of the nanotubes may elicit a signal. 
Hydroxypropyl cellulose, having a mass of 0.791 gm, was 
dissolved in a solvent of 25 ml chloroform for coating the 
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nano tubes to detect presence of HC1. In each case, an aliquot 
of 5 ml polymer solution was drop -deposited onto the 
SWCNT network in FIG. 1A to coat the corresponding 
SWCNTs. 

In one version of the system shown in FIG. 1, eight of the 
IDEFs (12A-1 through 12A-8) are doped with mono-layer 
protected clusters (“MPCs”) ofAu, eight of the IDEFs (12B-1 
through 12B-8) are coated with hydroxypropyl cellulose, and 
sixteen of the IDEFs (12C-1 through 12D-8) are “bare” 
CNTs. In another version, each IDEF, or subsets including 
more than one IDEF, is coated or doped with a (different) 
selected coating or dopant material so that the presence (or 
absence) of up to 31 or 32 different target materials can be 
tested. 

The electrical current through the sensors, at a constant 
voltage of 1 Volt, was monitored as different concentrations 
of chemicals, such as chlorine (Cl 2 ) and of hydrochloric acid 
(HC1) vapor were sequentially introduced to the sensor’s 
environment. A voltage difference of less than 1 Volt (or 
greater, if desired) can be used here. A computerized gas 
blending and dilution system, Environics 2040 (Environics, 
Inc. Tolland, Conn.), was used to create different concentra- 
tion streams with a steady flow of 400 cc/min during both 
exposure and purge periods. A gas cylinder containing 98.3 
ppm Cl 2 gas balanced with nitrogen, and a gas cylinder con- 
taining 478 ppm HC1 with nitrogen, were used as the source 
gases. Nitrogen was used both as the purge gas and as the 
balance gas for creating low concentration test samples. The 
test sample concentrations were 1, 2 and 5 ppm for Cl 2 gas 
and 5, 10 and 40 ppm for HC1 gas. The electrical signal 
(current) was collected using a semiconductor parameter ana- 
lyzer HP4155B (Agilent, Palo Alto, Calif.). Other equivalent 
electrical parameters, such as conductance or resistance, can 
be used as a response value. In trials involving heating, a 
thermal controller, Micro-Infinity ICN77000 Series Control- 
ler (Newport Electronics, Inc., Santa Ana, Calif.) with a ther- 
mocouple, maintained a constant temperature for the sensor 
operation. Additionally, a vacuum pump and an ultraviolet 
lamp of wavelength 254 nm were employed on occasion to 
accelerate the recovery of the sensors between tests; other 
ultraviolet wavelengths, such as 300 nm and 360 nm, can also 
be used to accelerate recovery. 

Carbon nanotubes do not sense presence of some gases and 
vapors due to the chemical and physical properties of CNTs 
as well as the nature of interaction between the gas molecules 
and nanotubes. SWCNTs have been found to detect presence 
of NH 3 and/or N0 2 , based on the charge transfer between 
these gases and SWCNTs. Our early tests indicated that pris- 
tine SWNTs do not respond at all when exposed to some 
industrial chemicals, such as chlorine and hydrogen chloride. 
It is important to get some observable response before one can 
do signal processing or pattern recognition for selective iden- 
tification. 

Carbon nanotubes coated with different polymers, such as 
chlorosulfonated polyethylene, hydroxypropyl cellulose, 
polystyrene, polyvinylalcohol, etc. used in commercial poly- 
mer based chemical sensors available for organic vapor detec- 
tion, can provide specific interactions with a chemical species 
of interest. As this chemical treatment aims to provide a 
specific interaction between the carbon nanotube matrix and 
specific gas molecules, the treatment can improve the selec- 
tivity while maintaining the high sensitivity expected of a 
nano sensor. 

Several polymer-coated carbon nanotube sensors have 
been tested for different toxic gases, such as chlorine and 
hydrogen chloride; for comparison, other gases, including 
oxides of nitrogen (NO p ; p=l-3), ammonia (NH 3 ), benzene, 
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nitrotoulene and acetone, have also been tested. FIG. 2 shows 
the room temperature sensor response for SWCNTs coated 
with chlorosulfonated polyethylene, for Cl 2 concentration 
pulses of 1 , 2 and 5 ppm as shown by the sample pulses. At 1 
5 ppm of Cl 2 , no response is observed, but significant response 
is seen at 2 ppm and above. The signal strength increases at 
higher concentrations. Such concentration-dependent signal 
response is a desirable attribute. By contrast, many state of the 
art sensors provide a flat signal regardless of the concentra- 
10 tion. The sensor results shown in FIG. 2 are consistently 
reproducible. However, the sensor recovery time is currently 
an issue, as is implicit in the times required («10 min) for 
reaction to change of concentration. From our previous work 
with NOp sensing, it is known that UV illumination of the 
15 SWNTs helps to drive out the adsorbents rapidly, thus accel- 
erating recovery when the source of the test sample is 
removed. 

Sensors with SWCNTs coated with hydroxypropyl cellu- 
lose have been tested for HC1 detection. FIG. 3 shows the 
20 sensor response, as a function of time, for HC1 concentration 
pulses of 5, 10 and 40 ppm. This sensor provides an observ- 
able response signal with 5 ppm HC1 present, with reasonable 
recovery times similar to those in FIG. 2. Again, the (equilib- 
rium) signal increases with the concentration of the analyte, 
25 and the sensor results are highly reproducible. Although the 
response and recovery times of this sensor for HC1 are 
improved relative to the results for chlorine, UV illumination 
and heating may be used here as well for further improved 
performance. 

30 A comparison experiment was conducted on sensors using 
pure, uncoated SWCNTs exposed to different gas and vapor 
analytes, with results shown in FIG. 4. The pure or uncoated 
SWCNT sensors showed no observable response when 
exposed to Cl 2 or to HC1, but displayed positive response 
35 signals, varying with concentration, for NO p , nitrotoulene, 
and benzene, and showed a negative response signal (oppo- 
site polarity) for NH 3 . Thus, the SWCNT sensor has some 
low level of discriminating power by itself to some, but not 
all, gases and vapors. It is clear from FIG. 4 that SWCNTs 
40 sensors are sensitive to NO^, nitrotoluene, benzene and NH 3 , 
as demonstrated previously. The opposite polarities of sensor 
response for NO^ and for NH 3 is believed to be due to the 
electron-acceptance and electron-donation characteristics of 
these gases. In relative terms, SWCNT sensors have lower 
45 sensitivity to benzene and to nitrotoluene vapors than to NO^ . 

A similar comparison experiment was carried out on the 
polymer-coated SWCNT sensors exposed to different ana- 
lytes, with results shown in FIG. 5. It is clear that addition of 
chlorosulfonated polyethylene coating to SWCNTs increases 
50 the sensitivity dramatically with a high level signal to 5 ppm 
Cl 2 compared to the pure SWCNT sensor, which shows no 
observable response even at 50 ppm of Cl 2 concentration. 
Although there is a signal for N0 2 from a chlorosulfonated 
polyethylene-coated SWCNT sensor, the sensing signal ratio 
55 of Cl 2 to N0 2 is much larger than 1 . Therefore, thi s sensor can 
be used in a sensor array to give a distinct response signal for 
Cl 2 (and, similarly, for HC1) for discrimination. 

These studies also show that both pure SWCNTs and chlo- 
rosulfonated polyethylene-coated SWCNTs do not respond 
60 to 100 ppm concentration of HC1 gas in nitrogen. Higher 
concentrations of HC1 were not tested as these high levels are 
not of interest for a nano sensor. In contrast, hydroxypropyl 
cellulose-coated SWCNTs respond to presence of HC1, but 
this sensor is also sensitive to N0 2 . Presence of the OH 
65 groups in the polymer coating may be responsible for the 
response signal differences in interaction with acidic gases. 
Because this sensor gives a significant response to HC1 that 
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other SWCNT sensors do not, the sensor can be used in a 
sensor array to provide a chemical signature that differenti- 
ates the HC1 gas from other chemicals. 

We have demonstrated a simple nano -chemical sensor 
using polymer-coated SWCNTs as the sensing medium. 
Because pristine (uncoated or unmodified) nanotubes do not 
respond observably to some chemicals of interest, it is impor- 
tant to explore coating or doping techniques to promote 
observable responses so that a broad application coverage can 
be ensured. We have found that the polymer coating enables 
selective sensing of chlorine and hydrochloric acid vapor at a 
sensitivity level of 5 ppm and above. It is important to recog- 
nize that coating or doping alone is unlikely to provide abso- 
lute discrimination. As with most sensors (of any size or 
exploiting any property change), pattern recognition tech- 
niques would be a valuable and necessary complement to 
provide discrimination. In that regard, the use of sensor arrays 
with multiple elements is an effective approach to chemical 
sensing, wherein the data from multiple sensors can be routed 
to a signal processing chip, integrated into the system, for data 
fusion and analysis. Advanced signal processing and pattern 
recognition techniques can be used to confirm (or refute) the 
assumed presence of a given species; in addition to the help 
from the selective coatings. Multiple sensing element arrays 
offer additional operational freedom when sensor recovery is 
slow and is a rate limiting process. Under such circumstances, 
a sensor would always be available while other sensors are in 
recovery mode. 

Using the results shown in FIG. 5, a procedure, shown as a 
flow chart in FIG. 6, can be implemented to detect presence of 
one, two or more target molecules, such as halogens, F 2 
and/or Cl 2 , or of one or more hydrogenated halogens, such as 
HF, HC1 or HBr, in an unknown gas G (which may contain 
N0 2 , NH 3 and one or more of the target gas molecules). In 
step 61, a known increment of a selected (possible) compo- 
nent of the gas G (reference component, such as N0 2 or NH 3 , 
or target component, such as Cl 2 or HC1) is added to and 
mixed with a portion of the gas Gto provide an augmented gas 
G\ In step 62, response values, V 0 (G;meas) andV 0 (G , ;meas), 
are measured for the uncoated SWCNT sensor In step 63, the 
SWCNT network is coated with a selected coating no. nl 
(e.g., chlorosulfonated polyethylene, hydroxypropyl, poly- 
styrene, polyvinylalcohol or another suitable (first) coating), 
is exposed to the gas G and (separately) to the gas G', and the 
corresponding response values V wl (G;meas) and V wl (G’; 
meas) are measured or otherwise provided. In step 64, the 
SWCNT sensor is coated (instead) with a different selected 
coating no. n2 (e.g., hydroxypropyl cellulose or another suit- 
able (second) coating), is exposed to the gas G and (sepa- 
rately) to the gas G', and the corresponding response values, 
V« 2 (G;meas) and V w2 (G';meas) are measured or otherwise 
provided. In step 65, the system uses a selected algorithm to 
estimate the concentration (e.g., in ppb) or concentration 
upper bound or concentration lower bound of one or more of 
the components believed to be present in the gas G (or in the 
gas G'). Two suitable algorithms for determination of gas 
component concentration are disclosed in flow charts in 
FIGS . 7 and 8, discussed in Appendices A and B, respectively. 

It is assumed initially in Appendix A that (i) the response 
value difference varies linearly with concentration difference 
of a single constituent that is present and (ii) the response 
value difference, in the presence of two or more gas constitu- 
ents in the gas G, is the sum of the response value differences 
of the single constituent gases. Linear response coefficients 
a y for the response value differences are assumed to be deter- 
mined experimentally or otherwise provided. As an example, 
assume that one reference gas (e.g., N0 2 or NH 3 ) plus first 
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and second target gases (e.g., Cl 2 and HC1), are suspected to 
be present in the gas G. Estimates of each of the concentration 
values c m0 for the initial (unaugmented) gas G are obtained 
from inversion of an MxM matrix equation relating these 
5 concentration values to response value differences for N coat- 
ings, where M (=2) is the number of gas components (refer- 
ence and target) believed to be present and N (^1) is the 
number of coatings (or dopings) used for the measurements. 

The approach discussed in Appendix A allows separate 
to weights, w„ and w'„ (or the same weight), to be assigned to the 
measurements of the initial gas and augmented gas. Prefer- 
ably, at least two of the weight values in Eq. (3) are positive 
(e.g., (w l5 w 2 ) or (w' 1; w' 2 ) or (w l5 w' 2 ) or (w' l5 w 2 )) for the 
example with N=2, and the relative sizes of the non-zero 
15 weights reflect the relative importance of the response mea- 
surements. If, as is likely, the four response measurements are 
believed to be equally important, one can choose 
w i =w 2 =w 'i =w, 2 = l* One can ignore one or two of the four 
measurements, in which event the corresponding weight val- 
20 ue(s) is set equal to 0. 

The response coefficients used in Eqs. (1 ) and (2) are not 

necessarily positive. For example, the response coefficient 
for the gas constituent N0 2 is positive for several of the 
SWCNT coatings used, while the response coefficient a 7; for 
25 NH 3 is observed to be negative for some of these coatings. 

Appendix B, and the corresponding flow chart in FIG. 8, 
discuss an approach for estimating an upper bound, or a lower 
bound, of a concentration value of a selected molecule (e.g., 
NO x or NH 3 or Cl 2 or HC1), again assuming that a response 
30 value difference varies linearly with a concentration value 
difference of the selected molecule. This bound is computed 
separately for each selected molecule. Appendices A and B 
apply to coated or doped CNT networks. 

Exposure of the coated SWCNT network to ultraviolet 
35 light can reduce the recovery time (normally ten hours or 
more) required to return the network to a substantially 
uncoated condition, by promoting accelerated detachment of 
the coating material from the SWCNT network. FIG. 9 
graphically illustrates measured response versus time with 
40 exposure to ultraviolet light for two different amounts of the 
same coating (sensors 1 and 2), indicating the improvement in 
recovery time where ultraviolet light is applied to lower the 
energy barrier to desorption of the adsorbed gas from the 
SWCNT network. The chemical and environmental time his- 
45 tories may cause the response curves to differ (FIG. 9). 

For some relatively small molecules, such as methane 
(CH 4 ), other hydrocarbons, and oxides of carbon (CO x ; x=l, 
2), an SWCNT network, doped with a transition element 
(“TE”,) such as Pd, Pt, Rh, Ir, Ru, Os and Au, can be used to 
50 detect presence of these molecules by detecting a change in an 
electrical parameter (conductance, resistance, current or volt- 
age difference) or response value associated with a path 
defined by an SWCNT network that extends between two 
electrodes having a controllable voltage difference or current. 
55 Some molecules, including nitrotoluene and phenol, are rela- 
tively strong electron donors and/or electron acceptors, and 
these molecules’ presence can be readily detected using 
“bare” or unmodified SWCNTs. Other molecules, including 
but not limited to methane, hydrocarbons and carbon oxides, 
60 manifest little or no electron donor or electron acceptor action 
so that monitoring an electric parameter value V of an 
unmodified or “bare” SWCNT network will, by itself, not 
indicate presence or absence of these molecules. 

Where CH 4 is adsorbed in a SWCNT/Pd matrix, the com- 
65 bination forms a weakly bound complex, such as Pd 6 (CH 4 ) -6 , 
where 6 is a relatively small positive number that need not be 
an integer. Methane, hydrocarbons and carbon oxides are 
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“greenhouse” gases and require detection capabilities in the 
ppb — ppm range to have much utility in environmental moni- 
toring. The sensing platform is similar to that illustrated in 
FIG. 1, where the SWCNT/TE compound serves as a current 
or voltage modifier whose electrical parameter value V 
changes when molecules of a selected target chemical, such 
as CH 4 or C m H w or CO x , are adsorbed on the SWCNT/TE 
compound, and the magnitude of the value V is monotoni- 
cally increasing (not necessarily linearly) with increase in the 
amount of the target chemical present. This approach for 
detecting presence of a taiget molecule is often at least one 
order of magnitude more sensitive at room temperature than 
is detection using catalytic beads or metal oxides, in part due 
to (1) nanoscale-induced charge transfer between TE-loaded 
nanotubes and CH 4 molecules and (2) a relatively large sur- 
face area per unit volume for SWCNTs, which permits 
enhanced adsorption of hydrocarbon molecules (e.g., CH 4 
and C m H„). 

Fabrication of a sensing platform for the SWCNT/TE net- 
work begins with sputter coating of about 1 0 nm thick Pd onto 
a pile of SWCNT powder. The TE-loaded or TE-doped 
SWCNTs are then dispersed into distilled, deionized water 
(e.g., 0.1 mg of SWCNT/Pd in 10 ml of the water). This 
solution is then sonicated and drop deposited onto interdigi- 
tated electrode fingers to create an electrical sensor with an 
initial resistance in a range of about 0.2-1 kilo-Ohm. Current 
through the network, with a 1 Volt difference, was monitored 
where 6, 15,30 and 1 00 ppm of CH 4 was present, using a gas 
stream flow of about 400 cc/min during exposure and during 
purge. Apart from preparation of the doped or loaded 
SWCNT network (as distinguished from coating a SWCNT 
network), the procedure for estimation of constituents present 
in a gas using a doped SWCNT network, is parallel to the 
procedure using a coated SWCNT network set forth in the 
FIG. 6 flow chart. 

Vacuum pumping and exposure of the network to ultravio- 
let light (5=254 nm) are used to reduce the recovery time of 
the SWCNT/Pd (i.e., removal of the CH 4 ) between tests (no 
recovery if these recovery acceleration procedures are not 
implemented). FIG. 9 graphically illustrates measured 
response versus time with ultraviolet light present, indicating 
the improvement in recovery time where ultraviolet light is 
applied to lower the energy barrier to desorption of CH 4 
molecules from the SWCNT/Pd network. 

FIG. 10 graphically illustrates sensor response (current in 
milli-Amps) versus time for each of the CH 4 concentrations 
of 6, 20, 60 and 100 ppm, and also graphically illustrates 
fitting of a logarithmic function, 

V/V^a In ( C) -6=ln{ C7 exp (6) } , (1) 

to relative response V/V 0 versus CH 4 concentration c, for 
different sensor networks. Resistance, as the electrical 
parameter, normally decreases with increasing concentration 
C of the constituent, while conductance, electrical current and 
voltage difference normally increase with increasing concen- 
tration c. The algorithm set forth in Appendix 1 and illustrated 
in Eqs. (A3)-(A6) can be applied to estimate concentration of 
one or more of the gases CH 4 , C m H„ and/or CO x , by replacing 
the concentration C 1 or C 2 or C 3 =C by the quantity 

x=ln{C7exp(6)}, (2) 

where the parameters a and b will vary with the particular gas 
constituent of interest. 

Methane, in the presence of the SWCNT network, may 
form a complex such as FI [Pd] *CH 3 . The H atoms in CH 4 tend 
to attract electrons from Pd, which in turn can obtain electrons 
from the SWCNTs to facilitate formation of the complex. 
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This behavior should also be manifest for some or all of 
similar transition metals, such as Pt, Ru, Rh, Ir, Os and Au. 
The detection lower limit for CH 4 at room temperature, using 
a Pd-doped SWCNT network, is estimated to be a few hun- 
5 dred ppb to a few ppm. This compares with a CH 4 detection 
lower limit Of 0.5-1 percent for conventional sensors, at 
temperatures T (min) i^450° C. 

FIGS. 11 and 12 graphically illustrate the environmental 
effects of varying temperature upon relative response V/V 0 
10 and upon response time, respectively, for SWCNTs plus Pd. 
The relative response in FIG. 11 appears to reach a relative 
minimum at a temperature T=T(min)«40° C. and to increase 
monotonically as IT -T(min)l increases. At room tempera- 
ture, T«20° C., the relative response is approximately three 
times the relative response at T=T(min). The response time 
(required for approximate equilibration of the response after 
exposure to a gas containing the target molecule) decreases 
monotonically with increasing temperature T and increases 
with the amount of load (coating or dopant) used to treat the 
20 SWCNT 

The SWCNT/TE sensor, with TE=Pd, has been tested at 1 5 
ppm and 30 ppm concentrations of CH 4 at gas temperatures of 
T=40° C., 80° C. and 150° C. The response parameter value 
(e.g., conductance or current) increases with increasing tem- 
25 perature, perhaps due to an enhanced catalytic effect of 
SWCNT/Pd binding with increasing temperature. 

Experimental results for other hydrocarbons, for CO x , for 
ketones and for aldehydes are qualitatively similar to those 
for CH 4 . 

30 ^ 

Effects of varying relative humidity (RH) on relative 
response, referenced to the response (V 0 ) at RH=0 percent, 
have been measured for several coating materials and doping 
materials, for a sequence of RH values. It is expected that 
relative response V/V 0 will decrease monotonically as the RH 
value increases, in part because the presence of a polar sub- 
stance such as water would interfere with, and partly mask, 
the change AW in an electrical parameter, for substantially all 
coating and dopant materials of interest. This expectation is 
borne out in measurements of V/V 0 , at T=T o =40° C., pre- 
sented graphically in FIG. 13, for the coating material 
hydroxypropyl cellulose, applied to a SWCNT, for RH=0, 15, 
30, 50, 70 and 90 percent. Each of these curves for V/V 0 can 
be approximated or estimated by a parametrized curve such as 

U/F 0 (RH;r 0 )=F e ((RH/RHo)" I ;a)«a-sech{(RH/ 

RH 0 r }+d-«) (3). 

where RH 0 is a reference RH value, m is a positive number 
and 0<a=il, each value being chosen for the particular coat- 
50 ing material or dopant material of interest. For sufficiently 
small values of the quantity RH/RH 0 , the parametrized curve 
in Eq. (3) is further approximated as 

55 V/V 0 (RH ;7o)« (4) 

a{l - (RH/RH 0 ) 2 "72} + (1 -a), = 1 - (a/2)(RH/RH 0 ) 2 A 

which is linear and decreasing in the quantity (RH/RH 0 ) 2m , 
60 and thus linear in the variable RH if 2m=l . More generally, 
the measured relative response value V/V 0 (RH;T 0 ) is sub- 
stantially monotonically decreasing in the value RH and 
resembles a trapezoid with a non-zero tail value, as illustrated 
in the generalized curve shown in FIG. 14. 

65 Introducing the dimensionless variable 

m=(RH/RH 0 ) wj , 


(5) 
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several interesting values of u can be identified on FIG. 15: (1) 
a knee value, u=w K1 , where the first w-derivative is changing 
most rapidly; (2) a maximum slope value, u=u MS , where the 
(negative) slope or u-derivative is a maximum value; and (2) 
a second knee value, u=u^ 2 (>u m ), where the first u-derivative 5 
is again changing most rapidly. The values are estimated from 
the identities 

(d 2 F e /du 2 ) wK =a ■ s ech 3 « { sinh 2 u- 1 }^=0, (6A) 

(a^^^^a-sechMcosh^-^KI, (6B) 10 

The value of the parameter a in the interval 0<a^i 1 is irrel- 
evant here. The solutions of Eqs. (6A) and (6B) are verified 

to be 

w^sinh - 1 { 1 )=ln{V 2 ± 1 } , (7A) 15 

ln{2±V3}, (7B) 

With appropriate choices of the parameter values RH 0 and m 
in Eq. (5), the experimentally observed locations of the values 
u MS and u^ 1 can be matched. 20 

An appropriate value of the parameter a can be determined 
as follows. Let u correspond to the value of u for which 
F e (u / )=f(0<f<l). From Eq. (3), this requires that 

a-sech(uj)+l-a=f (8A) 25 

sech(Wf)=l-(l-y)/a), (8B) 

which requires that a lie in a reduced range, l-f<a^l . Equa- 
tions (6A), (6B) and either (7 A) or (7B) can be used to 
estimate RH 0 , m and a. 30 

If the approximation in Eq. (3) is adopted, the zero point 
relative response V/V o (RH=0) for a particular coating or 
dopant material can be compensated for the presence of mois- 
ture (RH>0) by a compensation factor such as 

35 

Uo(RH=0)=U(RH>0)/{«-sech{(RH/RH o )" l }+(l-fl)}. (9) 

Other approximations, replacing the sech(u) function in Eq. 

(3) by sec{u}, or by exp{-u 2 /u 0 2 }, for example, can also be 
used here for compensation. 

For a given coating or doping material and fixed tempera- 40 
ture, the measured relative response V(RH;meas)/V 0 , as rela- 
tive humidity RH is increased over a sequence of values, can 
be compared with corresponding reference values V(RH;ref; 
h)/V 0 for each of a plurality of candidate gases (h=l, . . . , H) 
to determine if a particular candidate gas is present. Appendix 45 
C sets forth an analytical procedure for determining if a target 
gas is likely to be present, from a comparison of measured 
relative response values V/V 0 for a variable environmental 
parameter, such as relative humidity, temperature or pressure. 

FIGS. 15-1, 15-2 and 15-3 illustrate sequences of measure- 50 
ments of a response parameter W (electrical resistance) for 3 
IDEFs (12A-5, 12B-7, 12C-1) of the 32 coated or uncoated 
SWCNT sensors shown in FIG. 1, as a function of N0 2 gas 
concentration (2, 5 and 1 0 ppm), for varying relative humidity 
values of 0, 30, 50, 70 and 90 percent. The curves for RH=0 55 
percent and RH=30 percent are substantially coincident in 
most, but not all, situations, which suggests that for RH no 
greater than a certain threshold (e.g., RH^ r «30 percent), the 
response in change AW=(R/R0)-1 is substantially the same 
as for RH=0 percent. For most, but not all, of the IDEFs in 60 
FIG. 1, the parameter change AW monotonically decreases as 
RH increases. For some of the IDEFs, the measured AW 
values do not vary in a simple manner with RH value at higher 
N0 2 concentration values. 

FIGS. 16 A, 16B, 16C, 16D, 16E and 16F graphically 65 
present time variation of the parameter W=electrical resis- 
tance for a coated SWCNT in the presence of N0 2 with time 
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varying concentrations for relative humidity values RH=0, 
0.3, 0.5, 0.7, 0.9 and 1.0, respectively. In each of FIGS. 
16A-16F, the time varying N0 2 concentration c(N0 2 ) fol- 
lows a pattern set forth in Table 1 


TABLE 1 


Time Variation ofNO : 

> Concentration. 

0 t < 800 (sec) 

c(N0 2 ) = 0 

800 ^ t < 1500 

c(N0 2 ) = -0.5 (ppm) 

1500 <2400 

c(N0 2 ) = 0 

2400 ^ t < 3000 

c(N0 2 ) = 1 

3000 ^ t < 3900 

c(N0 2 ) = 0 

3900 ^ t < 4500 

c(N0 2 ) = 2 

4500 ^ t < 5300 

c(N0 2 ) = 0 

5300 ^ t < 6000 

c(N0 2 ) = 5 

6000 ^ t < 6900 

c(N0 2 ) = 0 

6900 ^ t < 7500 

c(N0 2 ) = 10 


In each of the graphs in FIGS. 16A-16F, the measured 
parameter value W=R/R0- 1 manifests a decrease whenever a 
(new) non-zero value of c(N0 2 ) is introduced, and the mag- 
nitude of this decrease increases monotonically with increas- 
ing time. However, for an intermediate time interval during 
which c(NO 2 )=0 ppm, the measured parameter value W 
undergoes a monotonic partial recovery toward a higher W 
value. This recovery, as measured by the slope of the value 
G(t) for an intermediate time interval during which c(NO 2 )=0 
ppm, is mildest for RH=0, and this recovery slope increases in 
value approximately monotonically as the RH value 
increases. 

For any one of the graphs in FIGS. 16A-16F, the temporal 
response of the W value in the different time intervals may be 
approximated as a sum of first order processes, 

W(r) = W 1 (rl )exp{- 1 (RH, c(N 0 2 ; l))(r-rl)} + W2(rl) ( 10 ) 

{1 - exp{-62(RH, c(N 0 2 ; l))(r - rl)}(c(N0 2 ; 1) > 0, t > rl) = 
W3(r2)exp{-63(RH)(r - r2)} + 

W4(r2){l - exp{-M(RH)(r - r2)}(c(N0 2 ) = 0, t > r2), 

W2(rl){l - exp{-£2(RH, c(N0 2 ; l))(r2 - rl)} = W3(r2), (1 1) 

where the N0 2 concentration is assumed to change from 
c(NO 2 )=0 at t<tl to c(NO 2 ;1)>0 at t=tl, and to change from 
c(NO 2 ;1)=0 to c(NO 2 )=0 at t=t2 (>tl). The exponent coeffi- 
cients bl, b2, b3 and b4 are monotonically increasing with 
increasing values of the relative humidity parameter RH. The 
coefficient value Wl(tl) is the initial value of the parameter 
W at t=tl , when the concentration changes from c(NO 2 )=0 to 
c(N0 2 )=c(N0 2 ;1); the coefficient W2(tl) is the asymptotic 
value of the parameter W, where the concentration c(N0 2 )=c 
(N0 2 ) is maintained indefinitely; and the coefficient W4(t2) 
is the asymptotic value of the parameter W where the con- 
centration c(NO 2 )=0 is maintained indefinitely. 

The measured relative response V/V 0 or parameter value W 
may also drift with elapsed time At, measured relative to a 
time at which the coated or doped SWCNT sensor was ini- 
tially prepared. If it is assumed here that the drift is a single 
first order process, the relative response, as a function of 
elapsed time At=t -t 0 , may be approximated as 

U(Ar)/U(g^-exp{(-ArAg}+(l-^) ? (12) 

where b is a non-zero value (positive or negative) and At 0 is a 
positive time increment. The quantity At 0 can be determined 
from the combination 
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V(2Ar)/V(r 0 ) - V(Ar)/V(r 0 ) = 


(13) 
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A? 2 =-A//{ln{ v{At)-b lexp(-Ar/Ar 0 i)+^l +62- 1}— In 
(62)}, 


( 28 ) 


exp{(- Ar / Ar 0 )}V (Ar) / V(r 0 ) - V(0)/V(r 0 ) 


The quantity b is then estimated from a relation 

M 1- ^Ary U(r 0 )}/{ 1 -exp{(-Ar/Ar 0 )}. (15) 

Drift with elapsed time At of the measured relative response 
value V/V 0 can then be compensated by computing a value 10 
V(At 2 ) in terms of a value at another time 

r(Af 2 )=F(A*i)+6 exp{ (- Af 2 /Af 0 )-exp { (- Af j/AZq) } , (16) 


U(r 0 )«U(AO/{^-exp{(-A?/Ar 0 )}+(l-^)}, (17) 

where the quantities b and At 0 are assumed known with rea- 
sonable accuracy. The correctness of the representation in Eq. 
(12) can be evaluated by forming 


R{ Ab=ln I V(2At)iV(t 0 )- V(At)IV(t Q )-\n\ V(At)/V(t 0 - V 

(0)/V(t o )\, (18) 20 

for two or more selected distinct positive values, At=Atl and 
At=At2, If the representation in Eq. (12) is substantially cor- 
rect, the ratio identity 


R(At l )/R(At 2 )=At l /At 2 (19) 2 5 

will be substantially satisfied for any distinct, positive values 
of At x and At 2 . 

* Where the temporal drift involves two or more different 
first order processes so that Eq. (19) is not substantially sat- 
isfied, for example, 30 

v(A t) •=V(At)IV(t 0 )>^hl exp{(-Ar/Ar 01 )}+62 exp {(-A?/ 

Ar 02 )}(l-M-^2), (20) 

the computations are somewhat more complex. For simplic- 
ity, it is assumed that At 02 »A/ 0 i and that Ati^At 01 so that 35 


exp { (- Ar/Ar 02 ) }»exp { (-At/At 0i ) } (21) 

and exp{(-At/At 02 )} can be ignored relative to exp {(-At/ 
At 0 i)}. This assumption should be verified a posteriori. One 


how forms the differences 40 

v(Af)-v(0)«6 1 (x : - l)+b2, (22A) 

v(2At)-v{At)^bl-x i (x i -l), (22B) 

v(3Ar)-v(2Ar)«^l-x 1 2 (x 1 -l) (22C) ^ 

x^exp^Ar/A^j, (22D) 

and the ratio 

50 

R' = |v(3Ar) - v(2Ar)}/{v(2Ar) - v(Ar)} (23) 

= {x\ - x\] I {x\ - x { }. 


Eq. (23) is a cubic polynomial in the unknown x l and can be 55 
factored into 

*i(*i-l)(*HD=o (24) 

for which the realistic solution is 

Xl =R, (25A) 60 

Ar 01 =Ar/ln{l/fl}. (25 B) 

From the preceding Eqs. (20A), (20B) and (23), one finds that 

b 1 ~{v(2At)-v{At)}IR \R - 1 ), (26) 65 

62={v(Ar)-v(0)}-{v(2Ar)-v(Ar)}/fl, (27) 


V(t 0 )=V(At)/{b 1 exp{(-Af/Af 01 )}+62 exp{(-A//Af 02 )} 

(1-61-62)}, (29) 

Recall that v(0), v(At), v(2At) and v(3 At) are measured quan- 
tities. 

FIG. 17 is a flow chart of a procedure for analyzing or 
compensating for drift with elapsed time At of a parameter 
value or relative response V/V(t 0 ). In step 171, the system 
receives at least three measured response values, V(t 0 ), V(t 0 + 
At) and V(t 0 +2At), for a response, where At>0. In step 172, the 
system determines if the response value varies with time, for 
example, by comparing two or more of the response values, 
measured at different times, At x and At 2 . If the answer to the 
query in step 172 is “no,” the system terminates the inquiry, in 
step 173, because the response value is substantially constant 
in time. 

If the answer to the query in step 172 is “yes,” the system 
determines if the response value V(t) varies according to a 
single first order process, in step 174, for example, using the 
analysis developed in connection with Eqs. (12), (18) and 
(1 9), or another suitable analysis. If the answer to the query in 
step 174 is “yes” so that the relative response set forth in Eq. 
(12) is substantially correct, the system determines the 
parameters b and At 0 , as developed in Eqs. (12)-(15), in step 
175. The system optionally compensates for temporal drift, 
for this single first order process, for example, as indicated in 
Eq. (17)), in step 176. 

If the answer to the query in step 174 is “no,” the system 
optionally assumes that Eq. (20) is substantially correct, in 
step 177, where At 02 «At 01 is assumed, and estimates the 
parameters bl, b2, At 01 and At 02 , for example, as discussed in 
connection with Eqs. (20) -(28), in step 178. The system 
optionally compensates for temporal drift for this multiple 
first order process, for example, as indicated in Eq. (29), in 
step 179. 

Appendix A. Estimation of Concentration of a Gas Com- 
ponent. 

FIG. 7 is a flow chart of a procedure for estimating con- 
centration of a gas component that is suspected to be present. 
This approach requires knowledge of certain response coef- 
ficients relating change of a response value to increase or 
decrease in concentration of a particular gas component G. In 
step 81, response coefficients a n m (n=l, . . . , N; m=l, 2, . . . , 
M) are estimated or otherwise provided for M single constitu- 
ent gases (numbered m=l , . . . , M; Mi^) (such as N0 2 , Cl 2 and 
HC1 or for any other three gases of interest), in the presence of 
SWCNTs individually coated with each of N selected coat- 
ings (numbered n=l, . . . , N; Ni^l), for example, with chlo- 
rosulfonated polyethylene (n=l) and (separately) with 
hydroxypropyl cellulose (n=2). 

For an SWCNT coated with coating number n, 

U„(G;meas)- U 0 (G;meas)=a„ 5l m c N0 2 +ct n x c cn + 

a n,3' C HCb (Al) 

where, for example, c N02 represents the N0 2 concentration 
(e.g., expressed in ppm or in ppb). The set of response coef- 
ficients {a n m \ m for different coatings, n=nl and for n=n2 
(*nl), will differ from each other, but each set is determined 
or estimated by measurement of the response value differ- 
ence, V w (G;meas)-V 0 (G;meas), of coated (n) versus 
uncoated (n=0), for each of the two (or, more generally, 
M^2) single constituent gases present in a known concentra- 
tion. For a single constituent gas N0 2 and no coating (n=0), 
for example, a 0 1 =0.034±0.002. 

In step 82, a known increment of one (or more) of the 
(suspected) constituent (e.g., N0 2 or Cl 2 or HC1), is added to 
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the gas unknown gas G to provide an augmented gas G\ In 
step 83, the response values, V 0 (G;meas) and V 0 (G';meas), 
for the uncoated SWCNT network (n=0), in the presence of 
the gases G and G’, are measured or otherwise provided. In 
step 84, the response values, V w (G;meas) and V„(G , ;meas), 5 
for the SWCNT network coated with the (single) coating 
number n, in the presence of the gases G and G', respectively, 
are measured or otherwise provided. 

In step 85, an error function €, defined by 

2e(x,y,z)=Z n w n -{ V n (G;meas)~ V 0 (G;meas)-a„ l c l -a „ 5 10 

c 2+a n 3C 3 } 2 +2. n w ■{ ^(G'jmeas)- UotG'jmeas)- 
<* n ^c\+&c x )-a na c 2 +a n jc 3 } 2 , (A2) 

is provided, where c l3 c 2 and c 3 refer to the concentrations of 
the reference molecule, the first gas molecule and the second 
gas molecule, Ac x is a known concentration increment of a 15 
selected one (Cj) of the reference molecule, the first gas 
molecule or the second gas molecule, added to the gas G to 
provide the gas G', and w M , and w'„ are selected non-negative 
weight values. The two sums in Eq. (A2) represent the con- 
tributions of the initial composition and the augmented com- 20 
position, respectively. These sums over n may include one, 
two, three or more coatings for which the response coeffi- 
cients are known. In this example, n=l, 2. 

The error function €(c ls c 2 , c 3 ) is to be minimized with 
respect to choices of the (unknown) concentration values c l5 25 
c 2 and c 3 . Differentiating € with respect to each of the vari- 
ables c : , c 2 and c 3 , in step 86, one obtains three coupled linear 
equations in these variables: 

30 

{(w! + W^aly + (w 2 + w ' 2 ) ol l } c l + ( A3 ) 

{(Wi 4- w\)ai,ia i?2 + (w 2 + ^2)02,102,2^2 + 

{(Wi 4- wi)0 U 0 U 4 - (vv 2 4 - W2)«2,102,3}C3 = 

{wi (Vi(G; meas) - V 0 (G; meas)) + 

Wy (V ! ( G 1 ; meas) - V 0 ( G ’ ; meas)) - 014 Ac?i )} 
fl i,i + {^2(V2(G; meas) - Vq(G; meas)) 4 - 

vi 4(V2(G / ; meas) - Vo(G'; meas)) - 02 ,iAci))0 2 ,i, 

{(W! 4- Wy )01,101,2 + (w 2 4- >^2)02,1 02,2 kl + ( A4 ) 

2/2 40 

{(w l + wi)0 u + (w 2 + W 2 )a 22 }c 2 + 

{0*1 4- wi)0 U 0!,2 + (w 2 + ^2)02,302,2)^3 = 

|wiVi(G; meas) 4- w'y(yy{G'\ meas) - 0i,iAci)}0i ; 2 + 

{w 2V2(G; meas) - Vo(G; meas)) 4 - 

W 2 {W 2 {G'\ meas) - Vq(G'; meas)) - 02 ,iAci)}02,2, 45 

10*1 + W^ai fl 0 1>3 + (w 2 + ^2)02, 102,3^1 + ( A5 ) 

{(wi + wi) 0 lf 2 0 i ,3 ++(w 2 + ^2)02,102,2)^2 + 
l(w 1 + >T)01,3 + (W2 + ^2)02, 3 ) c 3 = 

{wi V 1 (G; meas) - Vo(G; meas)) + 50 

(Vi (G'\ meas) - Vo(G'; meas)) - 01,1 Aci)} 

01,3 + {vv 2 V2(G r ; meas) - Vq(G; meas)) + 

W2CV2G'; meas) - 

Vo(G / ; meas)) - 02 ,iAci)}0 2 ,3- 

55 

In step 86, Eqs. (A1)-(A5) in the unknowns c l5 c 2 andc 3 are 
determined, using standard matrix inversion techniques, after 
verification that a 3x3 (more generally, MxM) coefficient 
matrix for the vector [c x c 2 c 3 ] tr has a non-zero determinant. 60 
These solutions, [c x c 2 c 3 ] tr , provide estimates of the concen- 
tration values of the corresponding chemicals in the gas G (or 
in the gas G’) in step 66 of FIG. 6. 

Preferably, at least two of the weight values in Eq. (A2) are 
positive (e.g., (w l5 w 2 ) or (w’ 1? w’ 2 ) or (w l5 w 2 ) or (w' 1; w 2 )), 65 
and the relative sizes of the non-zero weights reflect the 
relative importance of the response measurements. If, as is 
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likely, the four response measurements are believed to be 
equally important, one can choose w 1 =w 2 =w , 1 =w’ 2 =l. One 
can ignore one or two of the four measurements, in which 
event the corresponding weight value(s) is set equal to 0. 

The response coefficients a w m used in Eqs. (1) and (2) are 
not necessarily positive. For example, the response coeffi- 
cient a z> . for the gas constituent N0 2 is positive for several of 
the SWCNT coatings used, while the response coefficient a iy 
for NH 3 is observed to be negative for at least one of these 
coatings (FIG. 4). 

More generally, where M reference gas components (num- 
bered m=l, . . . , Ml) and target gas components (numbered 
m=Ml+l, . . . , M1+M2=M) with unknown concentrations 
are believed to be present and N coatings (numbered 
n=l, . . . , N), the error function € (analogous to Eq. (A2)) is 
defined by 

2 e(cy , . . . ,c M ))=2„w fJ -{F P fJ (G;meas)-F 0 (G;meas)- 
^ m a n ^c m } 2 +2. n w' n -{ F„(G';meas)- U 0 (G';meas)- 
2 ra 0„’i(c ra +AcJ} 2 , (A6) 

where one, or more then one, concentration value c m is aug- 
mented by a known amount Ac m . The error function € is 
minimized by differentiation with respect to each of the 
unknown concentration values c m . This yields M coupled 
equations 

2 w w w -{r„(G;meas)-F o (G;meas)-2 ra 0„^c OT o}0/ J ^o+ 

Z n w'A V„(G';meas)- V 0 (G , ;mesLs)-i: m a n ^(c m0 + 

^c mO )}a n ^ o =0, (A7) 

for index values m0=l, 2, . . . , M. These can be restated in a 
matrix format as 

2 M (w„-+w' M ){2 ra 0„, m c w o}0„, ra o} = 2„w„-{U„(G;meas)- 

Fo(G;meas)}+2„w'„-{2 w 0„^Ac wO )}0„^o, 

+2„w'„-{ r„(G';meas)- r 0 (G';meas)}. (A8) 

After verifying that the determinant of the MxM matrix of 
coefficients for the quantities c mo in Eq. (A8) is non-zero, this 
MxM matrix can be inverted to determine estimates for the 
concentration values c m0 (m0=l, . . . , M). These concentra- 
tion value estimates will depend, in part, upon the relative 
values chosen for the weight values w w and w' w for the coat- 
ings. Where one or more of the reference molecule concen- 
tration values c m0 (m0=l, . . . , Ml ) are known in advance, the 
estimates for these reference concentration values can be 
compared with the corresponding known values to evaluate 
the likely accuracy of the remaining estimated values. 

The approach set forth in this Appendix A can also be used 
to estimate an initial concentration value c m0 where the CNT 
network is doped or otherwise loaded, rather than being 
coated. 

Appendix B. Determination of Bound on Gas Component 
Concentration. 

A second algorithm does not require provision of a large 
number of response coefficients but only seeks to deter- 
mine if a particular target molecule is present in at least a 
selected concentration. For a selected coating, such as chlo- 
rosulfonated polyethylene or hydroxypropyl cellulose, on the 
CNT, a measurement of the response value difference AV=V 
(coated) -V(uncoated) is taken for modified gases, G'(l) and 
G'(2), where each of two distinct supplemental concentration 
values, A^mO) and A 2 c(m0), respectively, for a selected 
molecule no. mO (e.g., NO x or Cl 2 or HC1) is added to the 
original gas G. The concentration value c o (m0) of the selected 
molecule present in the original gas G is unknown, and the 
configuration of the CNT network is unknown. It is assumed 
that the response value difference AV increases approxi- 
mately linearly with the concentration difference Ac(m0) of 
the selected molecule so that 
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Ar 1 (mO)=v 0 +v 1 -(c 0 (mO)+A 1 c(mO)), 
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(Bl) 

AV 2 (mO)=v 0 +Vy-(c 0 (mO)+A 2 c(mO)). (B2) 

A molecule mO should be chosen for which IA 2 c(mO)-A 1 c 
(mO)l is at least equal to a selected positive threshold. The 
quantities v 0 , v x and c 0 (mO) are then related by the equations 

Vi = (AV 2 (mO) - AVi(m0))/{A 2 c(m0) - Aic(mO)}, (B3) 

vo + vi ■ co(mO) = AVi (mO) - vi ■ Ai c(m0) (B4) 

= AV2(mO) - vi ■ A2c(m0), 

and v 0 and \ l are determined, in part, by the CNT network 
configuration (assumed fixed and reusable) that is present. 
Where, as is likely, v o i^0, one infers that the initial concen- 
tration value c 0 (mO) for the molecule mO is limited by 

c o (m0)={AU 1 (m0)-v0}/v 1 -A 1 c(m0)^{AU 1 (m0)/v 1 }- 

A^mO) (B5-1) 

or 

c 0 (mO)={A V 2 (mO)~ v0}/v 1 -A 2 c(m0) ^ {A V 2 (mO)/v l }- 

A 2 c(m0) (B5-2) 

Equations (B5-1) and (B5-2) provide[[s]] an upper bound for 
the quantity c 0 (mO). Where it is known that the coefficient v 0 
is non-positive, Eqs._(B5-l) and (B5-2) can be inverted to 
provide lower bounds for the concentration: 

c 0 (mO)={A Vy (wO)-vO }/v 1 -A 1 c(mO) ^ (A Vy {mO)!vy }- 

A 1 c(m0) (B6-1) 

or 

c 0 (mO)={A V 2 (m0)-v0 }/v 1 -A 2 c(wO) §? (A V 2 (m0)/vy }- 

A 2 c(m0) (B6-2) 

This approach does not provide a direct estimate for the 
quantity C 0 (mO), only an indication of whether the molecule 
mO is or is not present in a concentration of no more than the 
right hand quantity in Eqs. (B5-1) or (B5-2). However, this 
approach does not require determination and use of the 
response coefficients a y . that are required for the putatively 
more accurate method set forth in Appendix A. The method of 
Appendix B can be used to estimate upper (or lower) bounds 
for concentration values C of one, two or more selected mol- 
ecules. 

FIG. 8 is a flow chart of a procedure for estimating an upper 
bound for a concentration value c 0 (mO) for a selected mol- 
ecule in the gas G. It is assumed that the response value 
difference AV varies approximately linearly with the concen- 
tration difference Ac(mO) of the selected molecule. In step 91, 
first and second (distinct) known increments, A^mO) and 
A 2 c(mO), of a selected molecule mO are added to a gas G to 
provide first and second augmented gases, G1 and G2. In step 


92, response value differences, 

A Vy(mO))=V(Gl ;meas)- U(G;meas), (B7) 

A U 2 (m0))=U(G2;meas)- U(G;meas), (B8) 

are measured or otherwise provided. In step 93, the coeffi- 
cient vl in an approximation for response value differences 

AU 1 (m0)=v o +v 1 -{c o (w0)+A 1 c(w0)}, (B9) 

A U 2 ( wO )= v 0 + v y{c 0 {mO) + A 2 c( mO ) } , (BIO) 

is determined according to 

Vy=(A V 2 (mO)-A Vy (mO))/ {A 2 c(m0)-Ayc(m0)}. (Bl 1) 


In step 94, the system queries whether the coefficient v 0 is 
likely non-negative. If the answer to the query in step 94 is 
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“yes,” the system estimates an upper bound for the initial 
concentration value c 0 (mO), in step 95: 

c 0 (mO)^{A Vy ( mO)/vy-Ayc(mO ) =A V 2 (m0)/vy-A 2 c(m0) 

5 (v o =0). (B12) 

If the answer to the query in step 94 is “no,” the system 
estimates a lower bound for the initial concentration value 
c 0 (mO), in step 96: 

10 c o (m0)^{AVy(m0)/vy-Ayc(m0) =AP 2 (w0)/v 1 -A 2 c(m0) 

(v o ^0). (B13) 

Appendix C. Effect of Varying Environmental Parameter. 

FIG. 18 is a flow chart of a procedure for estimating con- 
centration of a candidate gas component, suspected to be 
15 present in a target gas, where a response parameter V is 
associated with a controllably varied environmental param- 
eter (“EP”), and measurements of the response V(G;meas) are 
available for a sequence of reference values of EP, {EP^}^. 
The environmental parameter(s) EP may be one or more of 
20 the following: relative humidity RH, temperature T, gauge 
pressure GP, ambient chemicals present AC, and elapsed time 
At since the test sample was prepared. 

In step 181, measurements of the response V(EP p ;meas) 
(p=0, 1 , . . . , P;P^ 1 ) are provided, where p=0 may correspond 
25 to EP equal to a reference value, such as RH=0 or 0.1, for 
example. In step 182, a sequence of reference values V(EP fJ ; 
ref;h) (h=l , . . . , H:Hi^ 1) is provided, for each of H candidate 
gas components and for each of P reference values of the 
environmental parameter EP. In step 183, the system deter- 
30 mines an error value 


p ] (Cl) 

2_j Wp|V(£!P p ; meas) - V \EP p \ ref; h)\ q K/z = 1, ■ ■ ■ , //), 


where w p is a non-negative weight associated with the envi- 
ronmental parameter value EP^ and q is a selected positive 
number (e.g., q=l or 2). If desired, the approximation F e ((EP/ 
40 EP 0 ) m ;a) in Eq. (3) may be sub stituted for the reference quan- 
tity V (EP^ ;ref;h) . 

In steps 184 and 185, a minimum error value (or values) 

e(min)=miii{e(/z=l), . . . , e(h=H)} (C2) 

45 

is/are determined (e.g., €(min)=€(/z =h0)), and €(min) is com- 
pared with a selected threshold value €(min;thr). If 

e(min) =e (/z =/z 0) e (min ;thr) , (C3) 

50 the system interprets these conditions as indicating that the 
candidate gas component(s) corresponding to h=h0 is/are 
likely present in the target gas, in step 186. If the condition 
(C3) is not satisfied, the system interprets this result as indi- 
cating that none of the gas component(s) hO is likely to be 
55 present in the target gas, in step 187. 

What is claimed is: 

1. A method of sensing presence or absence of a specified 
gas molecule in a target gas, the method comprising: 

60 providing, in a chamber, a target gas that may contain up to 

H candidate gas molecules present with unknown con- 
centrations c h , numbered h=l, . . . , H (H^l), first and 
second interdigitated electrodes, connected to at least 
one of (i) a voltage source having a controllable voltage 
65 difference and (ii) a controllable current source, and a 
path Q that includes a plurality of carbon nanotubes 
(CNTs) and that extends between the first and second 
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interdigitated electrodes, where the path Q is loaded 
with at least one loading chemical and is exposed to the 
target gas; 

providing, within the chamber, an environment having a 
sequence of values EP^ (p=l , . . . , P;Pi= 1 ) of at least one 5 
environmental parameter, chosen from among tempera- 
ture, relative humidity, pH and gauge pressure; 
measuring, between the first and second interdigitated 
electrodes, in a presence of the target gas, an electrical 
parameter value (“response value”), denoted V(EP p ; to 
{c^} ;meas), of at least one of electrical current, electrical 
conductance, voltage difference and electrical resis- 
tance, where the target gas experiences an environment 
having at least one of the environmental parameter val- 
ues EP p ; 15 

providing reference response values V({EP };{c A (ref)}; 
ref) for a sequence of reference values of known con- 
centrations, {c^(ref)}={c A=1 (ref), . . . , c A=// (ref)}, for the 
path Q exposed to an environment having the environ- 
mental parameter value or values {EP^}; 20 

approximating the sequence of electrical parameter values 
V({EPp}; {c^(ref)};ref) as linear combinations, 

«F(EP^;{cJ;approx)=tf 0 +tfiCi+. . . +a H c H ,~V( E? p ; 

{c h ( ref)};ref) 

25 

of the concentration values c h of the candidate gas molecules, 
where a 0 , . . . , a H are known response coefficients that may 
depend upon one or more values EP^ of at least one of the 
environmental parameters that the candidate gas molecules 
experience, and a 0 represents a response coefficient when 30 
none of the candidate gas molecules is present; 

estimating an error value €({c^}) for the candidate gas 
molecule concentrations c h , the error value being 
expressed as 


H P 

£(k/J) = Yj XI w p\ y ti EP p}’ tad; approx) - V ({EP p }; {c h }\ meas)\ q , 

h = 1 p = 1 

40 

where w^ is a non-negative weight value that may vary with 
estimated relative importance of the environmental parameter 
values EPp and q is a positive number; 

determining or providing concentration values c^=c^(min) 
(h=l, . . . , H) in the approximate response values 45 
V({EPp}; {c^}; approx) that provide a minimum error 
value, € ({ c h } ;min)=€( { c^=c^(min) } ) , for the error value 

*({c*}); 

comparing the minimum error value €({c A };min) with a 
positive threshold error value €(min.thr); 50 

when e({c /? } ;min)^€(min;thr), interpreting this condition 
as indicating that at least one of the H candidate gas 
molecules, designated as h=hO, is present in the target 
gas, and that the concentration value of at least one 
candidate gas molecule that is present in the target gas is 55 
substantially equal to c A=A0 (min); and 
when e(c^;min)>6(min;thr) for all h=l, . . . , H, interpreting 
this condition as indicating that none of the H candidate 
gas molecules is likely to be present in the target gas. 
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2. The method of claim 1, further comprising 
choosing said exponent q equal to 2; 

partially differentiating said error value €({c^}) with 
respect to each of said concentration values c h 
(h=l, . . . , H) and setting each of the resulting partial 
derivative equations equal to 0, to obtain a system of 
non-homogeneous equations in said concentration val- 
ues c h ; 

obtaining solutions, c^=c A (sol), of the system of non-ho- 
mogeneous equations; and 

interpreting at least one of the solutions c^(sol) as said 
value c A=A0 (min) that would provide a minimum value 
for said error value €({c^}). 

3. The method of claim 1, further comprising choosing said 
weight values w^ to have substantially identical, positive val- 
ues. 

4. The method of claim 1, further comprising choosing said 
candidate gas molecule from a group consisting of NO x , NH 3 , 
acetone, benzene, nitrotoluene, F 2 , Cl 2 Br 2 , 1 2 , HF, HC1, HBr, 
HI, FC1, FBr, FI, CIBr, C1I, BrI, an alcohol, an aldehyde and 
a ketone. 

5. The method of claim 4, further comprising: 
drawing said loading chemical from the group of loading 

chemicals consisting of chlorosulfonated polyethylene, 
hydroxypropyl cellulose, polystyrene and polyvinylal- 
cohol; and 

providing said loading chemical by at least one of (i) coat- 
ing at least one of said plurality of said CNTs in said path 
Q with said loading chemical, and (ii) doping at least one 
of said plurality of said CNTs in said path Q. with said 
loading chemical. 

6. The method of claim 1, further comprising choosing said 
candidate gas molecule from a group consisting of at least one 
hydrocarbon and at least one oxide of carbon. 

7. The method of claim 6, further comprising: 
drawing said loading chemical from the group of loading 

chemicals consisting of Pd, Pt, Rh, Ir, Ru, Os and Au; 
and 

providing said loading chemical by at least one of (i) coat- 
ing said loading chemical on said at least one of said 
plurality of said CNTs in said path Q, and (ii) doping at 
least one of said plurality of said CNTs in said path Q. 
with said loading chemical. 

8. The method of claim 1, further comprising: 
exposing at least one of said plurality of said CNTs in said 

path Q to light having an ultraviolet component for a 
time interval having at least a specified time interval 
length At(min); and 

allowing at least one of said candidate gas molecules of 
said target gas to become desorbed from at least one of 
said plurality of said CNTs in said path Q. 

9. The method of claim 1, further comprising including on 
at least one of said first and second interdigitated electrodes a 
coating of at least one of Ti and Au, having a thickness in a 
specified range, deposited on an exposed surface of a sub- 
strate that is substantially nonconducting electrically. 

10. The method of claim 1, further comprising choosing 
said CNTs to include at least one single wall carbon nanotube. 





